Dielectric Anomalies of Solid CO and N 2 in the Audio Frequency Range 
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We report the first audio frequency dielectric constant measurements of CO and N2 in their solid 
phases. We have observed several new features, including (1) strong hysteresis effects above an 
onset temperature Th — 42 K in the /3-phase of pure N2, (2) absence of the expected short range 
antiferroelectric dipolar ordering in CO down to 4.2 K, and (3) anomalous temperature dependence 
of the dielectric constant in the a-phases of both CO and N2 . Quantum mechanical treatment of the 
molecular rotation explains some of the observed anomalies in the a-phase, but the strong hysteresis 
indicates that pure geometrical frustration plays a significant role in the /3-phase. 
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Audio frequency measurements of both the real and 
imaginary components of the dielectric susceptibility of 
solid CO and N2 have been carried out for the first time 
at temperatures ranging from their melting temperatures 
down to 4.2 K. These simple molecular solids have been 
a subject of intensive experimental and theoretical stud- 
ies during the last three decades ( |Q-||,|l@jl3,0 and the 
references therein). The interest in these systems derive 
from the fact that they serve as prototypes for studying 
the underlying physics of the simplest glass formers. The 
interactions in these solids are highly frustrated because 
the symmetries of the interactions and that of the lat- 
tice geometry are incompatible and as a result new phe- 
nomena occur. The molecules interact via electrostatic 
quadrupole moments which lead to a complex Pa3 anti- 
ferro-orientational ordering at low temperatures. The 
long range ordering is however fragile, and with the intro- 
duction of relatively small amounts of disorder (e.g. by 
replacing some N2 molecules with non interacting spher- 
ical Ar atoms) long range order is lost, and quadrupolar 
glass states with randomly frozen orientations are ob- 
served. These systems are believed to be analogous to 
the spin glasses and are particularly important because 
the interactions are short range and well known, and the 
order parameters and molecular dynamics can be studied 
directly. 

Dielectric constant/loss measurements in these sys- 
tems are particularly important because they directly 
probe the orientational ordering. However, only measure- 
ments at microwave frequencies have been performed in 
the case of solid N 2 [f| . The molecular rotations in these 
solids are strongly hindered, and as a consequence the 
relevant time scales for collective molecular reorientation 
are expected to be much larger than those for an ideal 
free rotor (~10 -12 sec). In fact, glassy behaviors have 
been observed for solid N 2 -Ar mixtures for which the 
characteristic reorientational time scale is ~10~ 4 sec [Q. 
Measurements in this more interesting audio frequency 
range have not been performed due to the lack of re- 
quired sensitivity. Nary et al. have carried out dielectric 
loss measurements of solid CO || in the audio frequency 



range, but to fully ascertain the nature of the dipolar re- 
orientations in the geometrically frustrated (3 as well as 
the ordered a-phases, one also needs to study the dielec- 
tric constant in the entire temperature range of the solid 
phase. Moreover, CO has an intrinsic dipolc moment; 
while this makes it more accessible for dielectric stud- 
ies because of the large change in its dielectric loss with 
temperature ||||, it also makes it more difficult to iso- 
late the much smaller contribution from the quadrupolar 
interactions which are responsible for the orientational 
ordering. We have developed a three terminal AC ca- 
pacitance bridge |t],|| with two parts per billion sensitiv- 
ity for measuring the real part of the dielectric constant 
H which allows us to study the collective orientational 
properties of solid quadrupolar molecular systems in the 
relevant audio frequency range. 

Figs. [I] and || show the real part of the dielectric con- 
stant e of CO and N 2 as a function of temperature T, in 
the 0.5 kHz to 20 kHz frequency range. The sensitivity 
and reliability of the apparatus can be judged from the 
observed jumps in e(T) at the structural phase transition 
at T a/3 , which is 61.57 K for CO and 35.61 K for N 2 . At 
this temperature the low T fee phase, the orientationally 
ordered Pa3 structure (or a phase) , undergoes a first or- 
der transition to a high T hep or f3 phase which does not 
support long range orientational order at any finite tem- 
perature. A second jump in e(T) occurs at the melting 
transition at T m , which is at 68.13 K for CO and 63.15 
K for N 2 . Note that while the figures show e up to three 
decimal place accuracy, our apparatus has a much higher 
sensitivity, so the error bars for e are negligible. However, 
as shown by the values of T a p and T m , the accuracy in 
the absolute value of temperature is within 0.3 K. 

In addition to the expected jumps, we observed several 
surprising features: 

1) For T < 25 K, s(T) for polar CO and non-polar 
N 2 are similar, including the dip in e(T) near 7.5 K. This 
shows the absence of any dipolar ordering in CO down to 
4.2 K. Also e(T) for CO did not depend on time (of the or- 
der of several days) which would otherwise indicate a slow 
thermal relaxation towards a dipolar ordered state. On 
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the other hand the dielectric response differs markedly 
from the conventional temperature dependence for a non- 
polar material described by the Clausius-Mossotti (C-M) 
equation. 
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FIG. 1. Dielectric constant of pure CO in the audio fre- 
quency range relative to e' , the value at 4.2 K. e' = 1.4211 at 
0.5 kHz. The inset shows the experimental and the theoretical 
results for the a phase below the dipolar freezing temperature 
Tp in greater detail. 

2) For T > 25 K and in the a phase, e(T) for CO is 
strikingly different from that of N2. The difference is due 
primarily to the dipolar contribution from CO. The sharp 
rise at Tp for CO can be attributed to the dipolar melting 
at Tp, above which dipoles begin to flip. We note that 
the results reported for N2-Ar-CO in || are qualitatively 
indistinguishable from our data for pure CO in this tem- 
perature range, except for a shift in Tp, where CO was 
used as a tracer to increase the dielectric response. This 
shows the importance of probing molecular quadrupolar 
systems without a dipolar tracer. 

3) The most surprising feature for N2 occurs in the /3 
phase, where it has a highly anomalous e(T). For ther- 
mal cycling below a temperature Th ~ 42 K, there is no 
hysteresis. In this regime, e(T) retraces a unique curve 
depending on the initial conditions, including the jump 
at T a p. However, if the system is taken above Th during 
a heating cycle (with the AC electric field present), the 



observed e(T) follows a different curve during cooling, for 
which s(T) is higher than the previous curve. The differ- 
ence between the two curves depends on how far above Th 
the system was heated. Fig. |] shows a typical hysteresis 
curve for a 1 kHz electric field of strength 1 kV/m where 
the system was heated to ~48 K along curve (1) and 
cooled along curve (2). When the system is heated to a 
temperature less than 48 K but above Th and cooled, the 
e(T) lies between curves (1) and (2). Once the system is 
cooled below Th, the corresponding e(T) curve remains 
reversible on thermal cycling, as long as the highest tem- 
perature remains below Th- Within 0.1 K accuracy, T a p 
remains the same for all curves. We should add that from 
these studies we cannot specify if Th is a sharp tempera- 
ture associated with the onset of hysteresis or a crossover 
region. We also note that the thermal resistivity data on 
pure N2 [I10| also has an anomaly close to T^. 
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FIG. 2. Dielectric constant of pure 14 N2 in the audio fre- 
quency range relative to e' , the value at 4.2 K on curve (1). 
e' — 1.4255 at 0.5 kHz. Curve (2) corresponds to Se/e' of the 
same sample heated above Th and cooled. 

4) If the system is left isolated (in the absence of an AC 
electric field) at a temperature above Th for a long time 
(of the order of several hours), e(T) retraces the lowest 
curve, independent of the initial conditions. 

5) No hysteresis was observed for CO in the entire a 
phase. The small window of temperature between T Q/ 3 
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and T m prevented us from carrying out similar measure- 
ments in the j3 phase of CO. 

6) The dielectric loss data for CO are similar to the 
results published elsewhere |||| but the loss data for N 2 
are featureless for all frequencies and temperatures stud- 
ied (not presented here). 

7) The frequency dependence of e at low temperatures 
can be collapsed on a unique curve when scaled by its 
value at 4.2 K. This is true for curves (1) and (2) and all 
other intermediate curves in Fig. ||. In general at 4.2 K, e 
increased with increasing frequency for both CO and N 2 
(~0.06% increase at 20 kHz from its value at 0.5 kHz), 
but we do not understand this frequency dependence. 
There is no frequency dependence for CO above Tp. 

While the phase of N 2 is novel, the behavior of e(T) 
in the a phase is also unconventional. In order to un- 
derstand the role of orientational ordering, we have de- 
veloped a phenomenological model based on the known 
temperature dependence of the orientational order pa- 
rameter for N 2 in the a phase. We find that although 
the molecules are large, we need to consider the discrete 
quantum mechanical orientational levels to explain the 
slow increase in e(T) in this regime. In contrast, the 
classical treatment gives only a linear temperature de- 
pendence. In the present model, the change in s(T) at 
T Q( g to be close to the experimental value requires 50% 
residual orientational ordering in the j3 phase. This is 
consistent with our conjecture that the ordering begins 
at Th > T a p. The same description holds also for CO up 
to the dipolar freezing temperature Tp. Above Tp the 
dipolar contribution overwhelms the orientational con- 
tribution. The sharp rise at Tp is associated with the 
onset of the flipping of the dipoles by 180°. However, 
the molecules are not yet free to rotate and are locked in 
the Pa3 structure. The subsequent decrease is therefore 
slower than the characteristic 1/T dependence expected 
for free dipoles. We find that a phenomenological model 
that accounts for the flipping without free rotation can 
describe the slower fall off qualitatively and predict the 
jump at Tp. For CO, the drop in e at T Q( g is consistent 
with the change in the order parameter with no resid- 
ual ordering in the (3 phase. Below we briefly sketch the 
phenomenological model for the a phases of CO and N 2 . 

The average potential energy of a diatomic molecule 
at a site in a crystal with Pa3 ordering can be written as 



V = CipP 2 (cos 9), 



(1) 



in which 9 is the angle between the molecular axis and 
the crystal axis (Z) at the lattice site and P 2 (cos0) = 
— ^ + f cos 2 (9) Q . The quantity ip is the average value 
of P 2 (cos#j) for the neighboring molecules j. C (~ —673 
K for CO and ~ —520 K for N 2 ) is a function of the 
molecular parameters and the lattice constants . The 
order parameter ip is known experimentally for both CO 
and N 2 in the a phase but little is known about the 



residual ordering in the (3 phase. If we write a = oti SO + 
Aa where a.i SO is the isotropic component and Aa is 
the anisotropic component of the electric polarizability 
of a single molecule in the condensed phase, the volume 
polarization can be written as [Ejfl: 



P = (l 



Xe 



)N(aiso + Aa((cos 2 (ft - 0)»)E e 



(2) 



where \e is the electric susceptibility, N is the num- 
ber density, £1 is the angle between the external electric 
field (Eext) and Z, and ((• ■ ■)) refers to both configura- 
tional and thermal average. In the present experiment 
at constant volume, N is a constant. Using the known 
values of the order parameter ip jllj], one can calculate 
e = l + Xe- However these calculated values of e increase 
linearly with T and do not agree with our experimental 
results. We note that when the separation of the adja- 
cent rotational energy levels of a rigid rotor are larger 
than kT, the classical Boltzman distribution used in the 
averaging in (2) may not be appropriate. In fact the sep- 
aration of the two lowest energy levels (AE/k) is ~ 93 K 
for N 2 and ~ 104 K for CO (setting ip = 1) @. From 
the Schrodinger equation for a rigid rotor with moment 
of inertia A and moving in the potential given by Eq. 1, 
the 9 dependence is given by 
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du> 



(1 



J 2 ) — 
dw 



Q 



fi — \oj 2 — m 2 (l — u 2 y 



= 0, (3) 



where u> = cos 9, fih 2 = ACip + 2AE, Xh 2 = 3ACip, and 
E is the energy. Eq. (3) is a spheroidal wave equation 
and the eigensolutions Q; m (A, ui) are oblate spheroidal 
wave functions. A is ~ —370 K for CO and ~ —275 K 
for N 2 For such large values of A, an expansion 

of the eigenfunctions in terms of Laguerre polynomials 
is appropriate. The eigenfunctions Qi m (X,ui) and the 
eigenvalues fj,i m (\,u>) for several of the lowest rotational 
states have been calculated using the method of Flammer 
Q . For various (ip, T) sets, ((cos 2 (f2 - 9))) is calculated 
first by determining Q; m and /i; m for several of the lowest 
eigenstates and the expectation values of (cos 2 (f2 — 9))i m 
for each of these states. Then 



((cos 2 (f7-0))) = i^ 7 ,(cos 2 (f7-0)), 



-Ej/kT 



(4) 



where Z is the partition function, "fj is the degeneracy 
of the jth state and Ej = (jij - ^§^j . The e cal- 
culated for N 2 using Eqs. (2,4) for various temperatures, 
is shown as the dashed curve in Fig. ^. This is in bet- 
ter qualitative agreement with the experiment. In order 
to obtain the correct jump at T a g, we need to set ip — 
0.5 in the ft phase. The same analysis has been car- 
ried out for CO with the known order parameter, and 
the resulting s(T) is shown as the dashed curve in Fig. 
□(inset). Again the agreement with the experiment is 
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qualitatively better. We mention two weaknesses of the 
above model. First, both N2 and CO show a small dip 
in e(T) near 7.5 K which remains unaccounted for in the 
present model. Second, if one takes the ai SO for a free 
gaseous molecule and number density N calcu- 

lated from the lattice constants ]Ti| , e' is about 15% and 
11% lower than the calculated value in the case of CO and 
14 N2 respectively. In Figs, [j] and || we have taken e' as the 
only adjustable unknown fitting parameter to compare 
the theory with the experiment. The lower experimental 
value compared to the theoretical estimate cannot be at- 
tributed to short-range antiferroelectric dipolar ordering 
in CO as suggested before || because N 2 also shows a 
similar behavior. However it can result from a lower po- 
larizability of a molecule in the condensed phase than in 
the gaseous phase, or from a lower number density for a 
powdered sample than for a single crystal. 

In the case of CO above Tp, the dipolar contribution 
to the dielectric constant should also be taken into con- 
sideration. Because the molecules are locked in the Pa3 
structure, the potential in the presence of an external 
electric field may be written as: 

V = CipP 2 (cos6) - ^ |E|/0)cos(r2 - 6). (5) 

The first part is the average potential energy of a single 
molecule as in Eq. (1). E is the internal electric field, fi 
is the intrinsic dipole moment, and / is some function of 
the order parameter. Qualitatively, when ip is unity (i.e. 
the molecules are completely ordered) the second term 
should be zero. When ip is zero (i.e. the molecules are 
free to rotate) it should simply be fi Q • E. For simplic- 
ity, we assume / to be of the form (1 — ip) T , consistent 
with the above two limits. The volume polarization can 
be obtained for the above potential as usual. Taking 
the same values for the isotropic polarizability and num- 
ber density which we used at 4.2 K, we find that for 
f(ip) = (1 — VO 1 / 8 the experimental data can be fit very 
well up to T a p (the dotted curve in Fig. 0). We note that 
the predicted change in e{T) at Tp is close to the change 
in the extrapolated value of the DC dielectric constant 
from experimental results within 0.8%. We have no ex- 
planation for the particular exponent in (1 — ^p) 1 ^ 8 , but 
the qualitative decrease which is slower than the charac- 
teristic 1 /T dependence of free dipoles can be understood 
within the model. 

In summary, we have made the first measurements of 
the dielectric constant of pure CO and N2 in the au- 
dio frequency range and observed several anomalous fea- 
tures. The low temperature behavior of both solids are 
very similar, but they are nevertheless very distinct from 
the conventional non-polar materials and can be qualita- 
tively understood in terms of a phenomenological quan- 
tum mechanical model. The strong hysteresis effects in 
the (3 phase of pure N2 show that the conventional pic- 
ture of the onset of ordering in these frustrated systems 



must be modified. Because the e(T) curve always re- 
traces the lowest curve in Fig. || whenever the system is 
left isolated at a temperature above Th for a long time 
and then cooled, we strongly believe that the above hys- 
teresis effects cannot be explained by the lattice defects 
such as dislocations or vacancies. However, the results 
are consistent with a departure from ergodicity due to 
trapping in a limited region of configuration space. As 
a consequence of the frustration of the interactions, the 
configuration space has a rugged landscape with many 
quasi-equal low energy minima separated by potential 
energy barriers. This phenomenon is common to a wide 
class of glass formers and other frustrated systems with 
the introduction of disorder. The most significant re- 
sult of this study is that the glassy effects are produced 
in a purely geometrically frustrated system without the 
introduction of disorder. This will be relevant to under- 
standing the co-operative behavior in other geometrically 
frustrated systems (pyrochlores, spinels, and nuclear an- 
tiferromagnets) . 

Stimulating discussions with Yasu Takano are acknowl- 
edged. We thank J. Xia, C. Parks, and P. Stachowiak for 
their help and gratefully acknowledge the support of all 
the personnel in our machine and electronics shops. This 
work is supported by a grant from the National Science 
Foundation No. DMR-962356. 



[1] T. A. Scott, Phys. Rep. 27, 89 (1976). 

[2] J. A. Hamida, N. S. Sullivan, and M. D. Evans, Phy. Rev. 

Lett. 73, 2720 (1994). 
[3] C. I. Nicholls, L. N. Yadon, D. G. Haase, and M. S. Con- 

radi, Phy. Rev. Lett. 59, 1317 (1987). 
[4] W. Kempinski and J. Stankowski, Low Temp. Phys. 21, 

74 (1995). 

[5] K. R. Nary, P. L. Kuhns, and M. S. Conradi, Phy. Rev. 

B 26, 3370 (1982). 
[6] Shang-Bin Liu and Mark S. Conradi, Solid State Com- 

mun. 49, 177 (1984). 
[7] E. D. Adams, Rev. Sci. Instrum. 64, 601 (1993). 
[8] M. C. Foote and A. C. Anderson, Rev. Sci. Instrum. 58, 

130 (1987). 

[9] S. Pilla, J. A. Hamida, and N. S. Sullivan, To be pub- 
lished. 

[10] L. A. Koloskova, I. N. Krupskii, V. G. Manzhelii, and B. 

Ya. Gorodilov, Sov. Phys. Solid State. 15, 1278 (1973). 
[11] V. G. Manzhelii and Y. A. Freiman, Physics of Cryo 

Crystals (AIP Press, Woodbury, NY, 1997). 
[12] Barbara C. Kohin, J. Chem. Phys. 105, 882 (1960). 
[13] B. A. Wallace, Jr., and H. Meyer, J. Low Temp. Phys. 

15, 297 (1974). 

[14] C. Flammer, Spheroidal Wave Functions (Stanford Univ. 

Press, Stanford, CA, 1957). 
[15] G. Maroulis, J. Phys. Chem. 100, 13466 (1996). 



4 



[16] K. M. Gough, M. M. Yacowar, R. H. Cleve, and J. R. 

Dwyer, Can. J. Chem. 74, 1139 (1996). 
[17] I. N. Krupskii, A. I. Prokhvatilov, and A. I. Erenburg, 

Sov. J. Low Temp. Phys. 1, 178 (1975). 



5 



